O ; Efficient formation of strongly bound ultracold 

^ ; cesium molecules by photoassociation with 

^ ; tunneling 

[ Mihaela Vatasescu^, Claude M Dion^ and Olivier Dulieu^ 

' ^ '1 ^ Institute of Space Sciences, MG-23, RO-77125, Magurele-Bucharest, Romania 

O ' ^ Department of Physics, Umea University, SE-901 87 Umea, Sweden 

^ Laboratoire Aime Cotton, CNRS, Batiment 505, Univ. Paris-Sud, 



C^ 



X 



91405 Orsay-Cedex, France 

E-mail: olivier . dulieuOlac . u-psud . f r 



C/3 

^^ ' Abstract. We calculate the rates of formation and detection of ultracold Cs2 

t^ . molecules obtained from the photoassociation of ultracold atoms through the double- 

^H ^l well Og (6S1/2 + 6P3/2) state. We concentrate on two features previously observed 

O ' experimentally and attributed to tunneling between the two wells [Vatasescu et al 

2000 Phys. Rev. A 61 044701]. We show that the molecules obtained are in strongly 

pg ' bound levels {v" = 5, 6) of the metastable a'^E^ (6S1/2 + 6S1/2) ground electronic state. 
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1. Introduction 

The control of elementary interactions between atoms or molecules in the gas phase is a 
long term concern for researchers, namely in order to find ways towards the full control 
of a chemical reaction, giving the ability to choose the reaction path unambiguously 
from a well-defined internal state towards a desired final state. In this context, laser 
cooling and trapping of atoms has opened an entirely new field of investigation, as they 
can be brought almost to rest, in a well-defined internal state. A spectacular example 
is provided by the observation of Bose-Einstein condensation in alkali gases 12 El , 
recently followed by the demonstration of quantum degeneracy in a fermionic alkali gas 

Slowing and cooling of molecules, despite the inherent difficulty caused by the 
complex internal molecular structure, represents an increasingly active research field 
with many achievements since the first observation of ultracold Cs2 molecules by 
photoassociation of laser-cooled cesium atoms [5^- Two approaches have been proven 
to be very efficient, namely Stark deceleration jB] and buffer gas cooling 7]. Other 
promising methods rely on phase space filtering of a molecular beam , on billiard-like 
collisions [Hj, or on a gas expansion out of a rotating nozzle ^01- A major breakthrough 
has been the observation of molecular Bose-Einstein condensates, using magnetic field 
tunability of Feshbach resonances [HI [T2]. Also, two groups recently measured for 
the first time rates for collisions between ultracold cesium atoms and cesium molecules 

dnun!. 

Among the various approaches, the photoassociation (PA) process is very attractive, 
as it starts from a pair of ultracold atoms which absorbs a photon to form an ultracold 
molecule in a short-lived excited electronic state [ISj. The main drawback is that the 
stabilization of the excited molecule into stable electronic states hardly provides an 
ensemble of molecules in a well-defined internal state, as it relies on spontaneous emission 
which populates many vibrational levels. Some attempts to use stimulated emission have 
been reported [T6t ITTj . but they were limited to the transition into high-lying levels. 
Indeed, due to the poor spatial overlap of the wave function of the photoassociated level 
— predominant at large interatomic distances — with the wave function of the lowest 
vibrational levels of the stabilized molecules, it is very difficult to create a significant 
proportion of ultracold molecules in their absolute ground state. A first indication of the 
possibihty to form ultracold molecules in their absolute ground state has been provided 
in jIHI for K2 molecules using a two-step photoassociation process. Very recently, a 
multistep excitation/deexcitation scheme has been set up to produce a fair amount of 
RbCs ultracold molecules in the v = level of the electronic ground state J19j . 

The purpose of this paper is to describe a possible way to fill these objectives, 
in the particular case of the formation of ultracold Cs2 molecules, using a single- 
step excitation scheme. It results from the cooperative action of photoassociation at 
large distances, tunneling and resonant spin-orbit coupling at short distances, and final 
stabilization by spontaneous decay towards a few low-lying vibrational levels of a stable 
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molecular state. We first invoked such a process for the interpretation of so-called "giant 
lines" observed in photoassociation spectra in [201, hereafter referred to as paper I. 
One of us also performed a time-dependent analysis of the tunneling process in such 
ultracold conditions, helpful for the determination of the characteristic times of the 
tunneling motion pi]. From the present calculations involving the available extensive 
photoassociative spectroscopy of Cs2, we demonstrate that the ultracold molecules are 
created with a very narrow distribution of vibrational levels, peaking at f " = 5, 6 of the 
theoretical potential curve of the Cs2 metastable triplet state. This could be the first 
example of the formation of ultracold molecules left mainly in a single deep vibrational 
level via single-photon PA. The possibility to use this mechanism to create the initial 
state of cesium dimers in ultracold collisions with cesium atoms has been discussed in 
ref.jIHl. 

We present our calculation of the photoassociation rate, of the cold molecule 
formation rate, and of the molecular ion signal resulting from the photoionization of 
the stabilized cold molecules. This emphasizes the role of the detection process in the 
interpretation of the strong intensity of these giant lines observed in the photoassociation 
spectrum. 

We start in section |21 by recalling the main facts concerning the tunneling effects in 
the photoassociation of Cs2. We then present the model used to calculate the molecule 
formation rates (section ^, followed by the results of the simulations for the Cs^ 
ionization signals (section |3]). Finally, concluding remarks are given in sectional 

2. Experimental and theoretical facts on giant lines and tunneling in the 
PA of cesium 

We briefly recall below the main results of paper I. The PA process between two ultracold 
cesium atoms in a magneto-optical trap (MOT), shown in figure ^ is written as 

Cs(6s2Si/2, F = 4) + Cs(6s2Si/2, F = 4) + h{u;B2 - Al) 

^ Cs2(0-(6s2Si/2 + 6p2p3/2; V, J)). (1) 

The pair of cold cesium atoms absorbs a photon detuned by Al to the red of the D2 
atomic transition frequency u;d2, to populate rovibrational levels of the external well of 
the double-well-shaped 0" molecular state correlated to the 6S1/2 + 6P3/2 dissociation 
limit (hereafter labeled Og (P3/2))- These levels are detected by two well-established 
methods: first through the fluorescence variations of the atomic trap induced by their 
spontaneous decay towards bound levels or towards the dissociation continuum of the 
lowest a^S+(6Si/2 + 6S1/2) state, second through the resonant two-photon ionization 
(R2PI) of the ultracold molecules formed in this metastable triplet state into CS2' 
molecular ions. The latter method delivers a very neat spectrum with excellent signal-to- 
noise ratio. Typical experimental conditions are temperatures between 20 and 200 /xK, 
a mean atomic density of about 4 x 10^° cm^'^ with a peak density of 10^^ cm~^, the 
number of atoms being estimated between 2 to 5 x 10^ atoms, and intensities of 50 
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Figure 1. Scheme of the photoassociation (PA) process and subsequent spontaneous 
emission (SE) of the 0" (P3/2) excited molecule to the metastable a'^E+(6Si/2 + 6S1/2) 
state. The 0~(P3/2) state is coupled to the Og (5D3/2) via spin-orbit interaction. 
Vibrational levels involved in the tunneling through the Oj (P3/2) barrier, as well as the 
deep a'^E^ vibrational levels which can be reached by spontaneous decay are pictured 
with lines (not to scale). The recorded experimental spectrum is recalled in the two 
insets, for both the Gi and G2 features, showing their large rotational structure, and 
the neighboring levels v^xt of the Og (P3/2) external well. 



up to 500 W/cm^ for the the photoassociation laser p2|- The full R2PI spectrum for 
the Og (P3/2) has been analyzed in jSHj, yielding a very precise determination of the 
corresponding potential curve through the Rydberg-Klein-Rees (RKR) analysis, and 
through an approach based on the asymptotic modeling of the atom-atom interactions 
^^. Vibrational levels from t>ext = to fgxt = 132 have been identified ^^ in the 
Og (P3/2) outer well, with a rotational structure (up to J = 4) clearly visible up to 
t;cxt = 72. 

In addition, two intense structures with a large rotational splitting (hereafter 
referred to as "giant hues" , following paper I, and labeled Gi and G2) are superimposed 
on lines associated with levels of the 0~ external well (figure Q). The most intense line 
within the Gi (G2) feature was assigned to J = 3 (J = 0, 1) at Al = —2.14 cm^^ 
(Al = —6.15 cm^^) superimposed on t>cxt = 103 (i^cxt = 80), while weak lines assigned 
up to J = 6 were detected providing a rotational constant B^^ = 137 ± 4 MHz 
{B^^ =243 ±8 MHz). 

With the help of a coupled-channel model, the Gi and G2 features have been 
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assigned in paper I to levels which tunnel through the potential barrier between the two 
wells of the 0~(P3/2) state, as schematized in figure [T] Their large rotational structure 
is induced by the vibrational motion in two coupled vibrational levels of the 0" (P3/2) 
internal well, which is coupled to several levels in the external well (only two of them 
are drawn in figure Q for better clarity). The tunneling effect, unusual for a heavy 
molecule like Cs2, is then an efficient mechanism to transfer the vibrational motion 
of the photoassociated state from large interatomic distances towards the inner zone: 
in contrast with the long-range molecular states usually reached by PA, spontaneous 
emission of the tunneling levels can stabilize the photoassociated molecules into low 
vibrational levels of the a^S^ metastable state, creating then ultracold molecules with 
a rather cold vibrational motion. Let us note that the internal state of the molecules 
created has not been probed experimentally yet. 

The 0g(P3/2) outer well is represented by the asymptotic model of [21]. As 
there is no available spectroscopic determination of the 0" (P3/2) inner well, this 
external well is matched to the curve computed by Spies and Meyer ^B], through a 
potential barrier whose height and position have been adjusted in paper I in order 
to reproduce the experimental observations. The potential barrier then culminates at 
2 cm~^ above the 6S1/2 + 6P3/2 dissociation limit. Furthermore, according to |26j . 
the Og (P3/2) state is coupled in the inner well region to the next 0" (6S1/2 + 5D3/2) 
state (hereafter labeled 0~(5D3/2)) through a non-adiabatic coupling generated by spin- 
orbit interaction. The 0" (P3/2) and 0" (5D3/2) curves then exhibit a well-locahzed 
avoided crossing around IOoq, which we transformed into a real crossing by linearizing 
the curves around the crossing to determine the corresponding standard Landau-Zener 
coupling parameters. We imposed a finite range for the coupling using the Gaussian 
form Aexp [— (r — roY/w'^] with tq = 10.02, w = 2, and A = 0.000258 (all in atomic 
units). We recall that the absence of resonant coupling between levels of the internal 
wells would result to a single tunneling level, instead of the two levels assigned in the 
experiment j20j . 

For each value of J, the radial Schrodinger equation is solved with the Mapped 
Fourier Grid Hamiltonian (MFGH) method ^27j , providing accurate vibrational energies 
Eyj, wave functions XvjiR), and rotational constants B^j = (x„j|^^/(2/ii?^)|x^j) (where 
/i = 121135.83 a.u. is the Cs2 reduced mass) of the two coupled state. It is well known 
that tunneling through a potential barrier is very sensitive to its shape and to the 
position in energy with respect to the top of the barrier. These effects are magnified 
here as the two tunneling levels of the Og (P3/2) and 0" (5D3/2) internal wells, hereafter 
referred to as Vqp and fsD, respectively, interact resonantly with the dense — but not 
continuous — energy level spectrum of the external well. Our calculations predict that 
the tunneling effect will be observable in the detection signal if the tunneling probability 
is maximum, which is achieved by adjusting the relative position of the two internal wells 
to create an almost half-and-half mixing of the f 6P and ^50 radial wave functions. Our 
model actually presents many tunable parameters for such an adjustment: indeed, the 
internal well of the potentials, as well as their coupling, are unknown experimentally, 
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so that their fine tuning is quite tedious. We also noted that the tunnehng effect is 
very sensitive to the value of J. In order to facilitate the convergence, we stopped 
the adjustment when we found maximal tunneling for an arbitrary J value. In other 
words, we also considered the rotational number J as a tunable parameter, which we 
will label as J in the following. An example of the resulting wave functions is shown 
in figure El, where the maximal tunneling effect is found for the value J = 5. This is 
clearly not in agreement with the rotational quantum number of the experiment, but we 
only need radial wave functions for the rate calculations of the next section, which will 
be evaluated without any J dependence. The corresponding numerical data is available 
on request from the authors. 

Even if in the adjustment procedure above we used J as a parameter, figure El 
shows that the selectivity of the tunneling process is very strongly dependent on the 
rotational level. The value of B^j defined above is very sensitive to the repartition of 
the wave function inside the two wells, and panel (a) in figure El refiects the pattern 
shown in figure El where the value for J = 5 is intermediate between the rotational 
constant of levels in the internal well (~ 0.008 cm~^) and of levels of the external 
well (~ 0.0005 cm^^). Similarly, the J = 3 line for Gi also refiects the extension of 
the wave function over both wells. The J = 5 line displays the same behavior, which 
can be understood since the rotational structure of Gi has the same magnitude as the 
energy spacing between vibrational levels of the external well, so that it is resonant 
with fext = 104 in our calculations. This aspect was not present in the model we set up 
in paper I and confirms what is observed experimentally. Indeed, due to the ultracold 
temperature, the s-wave regime is expected to dominate the initial collision between the 
two Cs atoms (yielding J = 0,1,2) with probably a small contribution of the p wave 
(yielding J = 1,2,3). The presence of higher values of J at such low temperatures is 
not yet fully understood |2S]. The relative intensities of the J = 0,1,2 experimental 
lines approximately refiects their 2 J + 1 degeneracy, and the J = 2 line has indeed the 
strongest intensity for most of the PA lines assigned to the levels of the 0~(P3/2) external 
well. The lines with larger J values are weaker, due to the small contribution of higher 
partial waves involved in the collision. We see that this hierarchy is not preserved for 
Gi and G2: the intense lines are associated with those rotational levels which are indeed 
resonant with a level of the internal well, i.e., J = 3 and J = 0, 1, respectively. 

3. Photoassociation and cold molecule formation rates for tunneling levels 



The metastable a S+(6Si/2 + 6S1/2) state is represented by the potential curve of 
matched beyond 20ao to the standard asymptotic expansion I] C„/-R" (n = 6,8,10), 
with Cn coefiicients taken from [SH]- The repulsive wall of the potential is slightly 
changed in order to reproduce a large and positive scattering length (SU]. We work 
here with ot = 2370ao, which satisfactorily reproduces the intensity envelope of the 
PA spectrum of j25, . The initial scattering radial wave function is computed through a 
standard Numerov integration at an energy E/k-sT = 150 /iK above the 6S1/2 + 6S1/2 
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Figure 2. Two-channel vibrational wave functions Xv .j{R) (panels b-e) associated 
with the radial motion in the 0" (P3/2) double well (in black) and in the 0" (5D3/2) 
well (in red). The corresponding potentials are drawn with full black and dashed red 
lines, respectively, to make more visible the extension of the wave functions within the 
corresponding potentials. The J — 5 value (panels b and c) allows maximal tunneling 
for the G2 resonance at —5.87 cm~^, almost degenerate with a level of the external well 
at —5.89 cm~^: the wave functions have significant amplitude over the barrier range. 
The J = case (panels d and e) is also shown for comparison: the energy matching is 
less favorable and the wave functions hardly penetrate the potential barrier. Panel (a) 
displays a typical radial wave function of the initial coUisional state for a temperature 
of 150 /iK, in the a'^S+(6Si/2 -I- 6S1/2) state, designed here with a scattering length of 
about 2370ao. 



asymptote, ignoring the hyperfine structure (figure |2t') • Let us note that the accuracy 
of the ax and C„ values is not crucial for the following calculation of the rates, which 
can be determined experimentally typically within a factor of 2. 

The photoassociation and cold molecule formation rates are evaluated according to 
the perturbative model already presented in previous papers |22| IHTI |32] • Briefly, the 
photoassociation rate Rpa{v, J; T) per atom (expressed in s~^) for an initial continuum 
state with energy E = ksT and wave function iPe{R) into a vibrational level v with 
a radial wave function Xj;j(0~;i?) of the coupled 0" states (where v stands for the 
quantum number v of the 0~(P3/2) external well levels, or for Gi or G2), at a detuning 
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Figure 3. Rotational constants Byj — {x^ j|?i (2/ii?^)|x5 j) for rovibrational levels 
of tiie coupled Og (6S1/2 + 6P3/2) and 0" (6S1/2 + 5D3/2) potentials, as function of their 
binding energy, (a) G2 structure; (b) Gi structure. The levels with an important 
tunneling effect are marked by arrows. 



A^ below 6S1/2 + 6P3/2 limit is expressed as 

/ 3 \ '^^'^ h E 



(2) 



where n^t is the atomic density and Ath = hJl/{3jj,kBT) is the thermal de Broglie 
wavelength. The atomic Rabi frequency 2K is related to the intensity / of the laser 
through K'^ =87"' ^^ere T/27i = 5.22 MHz is the natural width of the 6P3/2 atomic 
level. At the PA wavelength considered, ApA, the saturation intensity Jq = nhcT/^SXppJ 
is 1.1 mW/cm . The angular factor A = 125/3888 includes hyperfine degeneracies, 
assuming an initial {F = 4) + {F = 4) hyperfine state. As mentioned above, we actually 
solve the Schrodinger equation for every chosen value of the rotational quantum number 
J taken as an effective parameter, so we don't include the degeneracy factor 2 J + 1 in 
the rate formulas. 

The two components of the PA level wave function Xci,j(Og ; R)) (see figure Ej) are 
denoted by x^j{R) and x^^(-R), such that 

OiE, V, J) ^ \{^Pe |x.,j(Og )) f = \{^E \xlj) + {i'E \xfj) f • (3) 

Following [HB], the rate Rmo\{^h) for cold molecule formation after PA in the v level is 
obtained by multiplying the PA rate with the branching ratio Rhri^-, J) of the v level 
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towards the bound levels f "(a^S+)) of the metastable triplet state, and neglecting here 
the R dependence of the dipole transition function for the spontaneous decay step, 

i?^ol(AL) = RpAi^L)Rhriv, J), (4) 

with 

RUv, J) = E |(x.,j(Og ) 0."(a^S+)) ' . (5) 



Figure 0^ displays the overlap between the continuum wave function and the bound 
states of coupled potentials, equation 0, lying in the energy interval between —9 and 
—1.5 cm"^, for both J = and J = 5. As expected, Gi and G2 have smaller overlap with 
the initial continuum than the external well levels, due to the weaker probability of wave 
functions to be localized at large distances. In contrast with the PA rate, the branching 
ratio of the Gi and G2 wave functions is enhanced due to their good localization at 
short distances (figure Eb). As it was often emphasized, a good production of cold 
molecules requires a favorable ratio between free-bound and bound-bound transitions. 
The product 0{E,v, J) x -Rbr(^5 J) is represented in figure Et. If for Gi the efficiency 
seems certain, we see that for G2 this balance is fragile and depends indeed on the 
J value which makes tunneling effective. On the right axis of this figure is reported 
the rate for cold molecule formation, equation Q, for typical experimental conditions 
nat = 10^^ cm~^ and a PA laser intensity / = 100 W/cm . It is found two times and 
eight times larger than the rate for the Vext levels for G2 and Gi, respectively. Even if in 
the experiment the intensities of the G2 and Gi lines are found to be comparable (see 
insets in figure Q), our results confirm that the cold molecule formation rate is indeed 
larger than the rate for the fext levels. The relative intensity of the G2 and Gi lines may 
be infiuenced by the two-photon ionization used for the detection (see next section). 
Let us note also that the rate computed for the fext levels is in agreement with the one 
reported in [22], where Rmoi ~ 0.2 s~^ was measured at a detuning of 7 cm~^, for a PA 
laser intensity / = 55 W/cm . 

The rates for formation of cold molecules in individual levels v" of the metastable 
a zj state, 



i?^oi(AL,i;")=^PA(AL) (x.,j(0-) 0."(a3S+)) , (6) 

are represented in figure El for the same range of detunings and same experimental 
conditions as in figure ^, for both the Gi and G2 features, considering the levels 
indicated by arrows in figure |^ corresponding to J = 5. Both vibrational distributions 
peak markedly at v" = 5, 6, confirming the efficiency of the PA into the tunneling 
resonances to produce cold molecules in very low vibrational levels. 

4. Modeling the ionization signal 

In cesium PA experiments ED], the cold molecules are detected via a two-photon 
ionization process, where the first photon is resonant with rovibrational levels of 
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Figure 4. (a) Overlap 0{E, v, J) between the continuum wave function and 
vibrational wave functions of coupled 0~ potentials lying in the detuning interval 
between —9 and —1.5 cm^^. (b) Branching ratio Rhr{v, J) of the same levels to the 
a'^E+ vibrational states, (c) Product of the two previous quantities (left axis) and 
corresponding cold molecule formation rate R^oii'^h) (right axis) assuming an atomic 
density riat = 10^^ cm^"^ and a PA laser intensity / = 100 W/cm . Note that in this 
panel the contribution of the nearby levels at —5.87 cm~^ and —5.89 cm^^ shown in 
figure 12 are added together into a single rate value around —5.9 cm~^. Both J — 
(black crosses) and J = 5 (red open circles) are shown. The levels corresponding to 
Gi and G2 are indicated. The arrows identify the levels for which we calculated the 
vibrational distribution of the produced cold molecules displayed in the figure |3 



molecular states correlated to the 6S1/2 + 5D3/2 dissociation limit. Therefore, the 
ionization signal is very sensitive to the wavelength chosen for the ionization laser. 
Following our previous studies of this process jHZl IHBI, we simulate the contributions 
of Gi and G2 to the CS2' ions signal obtained after PA in a 0~(6Si/2 + 6P3/2) level 
detuned by Al, spontaneous emission towards f"(a^S+) bound levels, and ionization 
with a photon of frequency z/jon of these levels through absorption into levels v' of the 
(2)^ng(6Si/2 + 5D5/2) potential. In this model, the second step of the ionization process 
is described as a uniform ionization probability of the v' levels into Cs^ ions, and will 
not be considered here. The ionization signal Sior^^Ai^) is expressed as 

2 



^ion(AL)=iVpA(AL)E^ionK') {X,,j{%) <PA^'^t 



(7) 
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Figure 5. Rate of formation of cold molecules in the vibrational levels w"(a^S+), via 
photoassociation into (a) the Gi and (b) the G2 states. In the latter case, rates for 
both levels contributing to G2 are given, i.e., the level at —5.87 cni~^ tunneling from 
the inner wells (black closed circles) and the external level tunneling from the outer 
well at —5.89 cm~^ (red open circles). The calculation is made for riat = 10^^ cm~^ 
and an intensity / — 100 W/cm for the PA laser. 



As in the cold molecules rate calculation, we neglect the R dependence of the dipole 
transition moment involved in the spontaneous decay step. In equation ((7j), the number 
of molecules A^pa(Al) accumulated during the photoassociation step is defined by 
multiplying the PA rate /2pa(Al) with the number of atoms ripA in the photoassociation 
area with the residence time tpA of the cold molecules within the trapping region: 

ArpA(AL) = i?PA(AL)npAtpA. (8) 

In the Cesium PA experiment, typical values for these parameters are ripA = 5 x 10^ 
and tpA ~ 10 ms. The excitation probability Pion{v") of the v" levels into vibrational 
levels v' of the (2)^ng(6Si/2 + 5D5/2) potential is calculated as 
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Figure 6. Cs^ ions signal (in arbitrary units) as a function of the frequency of the 
ionization laser, calculated with equation for detunings Al corresponding to Gi 
(right column) and G2 (left column). Signals corresponding to photoassociation in 
vibrational levels located close to Gi and G2 and belonging to the external well of 
the Og (6S1/2 + 6P3/2) potential are also shown. The energy of the 6S1/2 + 5D5/2 
dissociation limit is taken as being 14597.08 cm^^ above the 6S1/2 + 6S1/2 origin. 



where D is the i?-dependent dipole moment for the a^S+ 
calculated in jSE], and 



f{^v"v' 



exp 



ln(2) 



{Ujj/iyi 



(Su) 



^ng(6s + 5D5/2) transition 



(10) 



is a parametric function accounting for the resolution of the experiment, estimated at 
Su ~ 30 GHz IH2], and u^ny/ is the transition frequency between v" and v' levels. 

The ion signals as a function of the ionization frequency z/ion obtained for Gi and 
G2 are displayed in figure El together with those corresponding to PA in neighboring 
vibrational levels fext of the Og (6S1/2 + 6P3/2) external well. As expected, large 
amounts of molecular ions are detected only for specific frequencies corresponding to 
the resonance condition in the first step. The signals coming from Gi and G2 are clearly 
more intense than those of the surrounding levels, in agreement with the experimental 
spectrum [20] , because of the increased efficiency for forming cold molecules in bound 
levels of the triplet state. The signal corresponding to Gi, having a maximum of 
0.03 in figure |B^, is five times bigger that the simulated G2 signal, whose maximum 
reaches 0.006 if we take the sum of the two contributing levels (figures |BK,b). This ratio 
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Figure 7. Same as figure we superimposed predicted ion signals from panels a, b, 
and e of this figure, slightly shifted from each other for better visibility. Note that the 
signal for Gi is shown at 1/lOth of its actual value. 



between the intensities of the simulated "giant lines" is larger than in the experimental 
spectrum where the ratio is only two I^Sj- This could be due to the fact that the 
ionization spectra for the tunneling resonances do not match perfectly: the differences 
in vibrational distributions w"(a^S+) of Gi and G2 (see figure EI) make it possible to 
chose an ionization frequency (e.g., huion ~ 14101 or 14117 cm~^) such that the ratio 
between the two ion signals is greatly modified, as shown in figure H 



5. Conclusion 

As demonstrated through several studies j31 1201 IHB] , cesium atoms are well-suited for 
ultracold molecule formation through photoassociation, mainly due to the peculiarity of 
some of the excited states of the Cs2 molecule. Of particular interest is the 0~ (S1/2+P3/2) 
state, which is composed of two potential wells separated by a barrier whose height 
culminates at an energy close to the dissociation asymptote. This allows for tunneling 
between the two wells, enabling a pair of atoms initially far apart to come close together 
in a photoassociated molecule. 

We have shown here that this tunneling process leads to strongly bound ultracold 
cesium molecules in the metastable a^S+(Si/2 + S1/2) ground electronic state. Indeed, 
it is mostly vibrational levels around v" = 5,6 that are populated, in contrast to what 
can be reached for non-tunneling states of the external 0" (P3/2) well {v" > 23 |22j). 
The tunneling efficiency is seen to be markedly dependent upon the rotational state J, 
in agreement with what was observed experimentally |2Uj . 

Finally, we have investigated that the experimental detection scheme, resonant two- 
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photon ionization, influences the observation of ultracold Cs2 molecules. Evidence would 
suggest that the detection process affects the relative intensities seen in the signal, as 
was the case previously |HS1- In other words, the molecules produced are not all ionized 
with the same efficiency, and it is therefore important to consider this step in the analysis 
of photoassociation experiments. 
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